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Abstract. The possible solitary wave solutions for a general Boussinesq (GBQ) type 
fluid model are studied analytically. After proving the non-Painleve integrability of 
the model, the first type of exact explicit travelling solitary wave with a special 
velocity selection is found by the truncated Painleve expansion. The general solitary 
waves with different travelling velocities can be studied by casting the problems to 
the Newtonian quasi-particles moving in some proper one dimensional potential fields. 
For some special velocity selections, the solitary waves possess different shapes, say, 
the left moving solitary waves may possess different shapes and/or amplitudes with 
those of the right moving solitons. For some other velocities, the solitary waves are 
completely prohibited. There are three types of GBQ systems (GBQSs) according to 
the different selections of the model parameters. For the first type of GBQS, both 
the faster moving and lower moving solitary waves allowed but the solitary waves 
with "middle" velocities are prohibit. For the second type of GBQS all the slower 
moving solitary waves are completely prohibit while for the third type of GBQS only 
the slower moving solitary waves are allowed. 
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1. Introduction 

The study of the Korteweg de-Vries (KdV) equation 

u t + Quu x + u xxx = (1) 

has been an interesting issue since the discovery of soliton. Its exact solution [|TJ can 

be used to describe the overtaking collision of soliton on a uniform layer of water, but 

the solution is only physically meaningful for the unidirectional soliton. All the left 

moving solitons with zero boundary conditions are prohibited. 

In Ref.||, three sets of model equations are derived for modelling nonlinear and 

dispersive long gravity waves travelling in two horizontal directions on shallow waters 
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of uniform depth. A good understanding of all solutions of these models are helpful 
for coastal and civil engineers to apply the nonlinear water wave model in a harbor 
and coastal design. Therefore, finding more types of solutions of these equations are 
fundamental interest in fluid dynamics. For the various different models, a point of 
central interest is to examine and compare the solitary wave solutions. For the case 
of plane motion, the three set of models can be rewritten as the Boussinesq class after 
omitting the higher order terms: 

(A) . The {u, v} system, - the depth- mean velocity basis 

1 

u t + uu x + v x = -u xxt ^ 
v t + [(1 + v)u\ x = 0. 

(B) . The {u, v} system, - the bottom variable basis 

1 

u t + uu x + v x = -u xxt 



v t + [(1 + v)u} x = -u xxx . 



1 

6" 

(C). The {it, v} system, - the surface variable basis 



(3) 



u t + uu x + v x = 

1 (4) 
vt+[(l + v)u] x = --u xxx . 

In the Boussinesq types of systems (HI)-®, the field v is wave elevation and u, u and 
u are depth-mean, bottom and surface velocities respectively. 

To study the systems (H)-(E1) uniformly, we can extended them to the following 
generalized Boussinesq system (GBQS): 

U t + UU X + V X = CiU xxt 

(5) 

v t + [(1 + v)u} x = c 2 u xxx . 

The system (|J) has been shown to be integrable and has Hamiltonian structures |J. 
The exact solitary wave and periodic wave solutions of (|4]) have also been obtained 
by many authors |5|]. Especially, Zhang and Li|| presented a theory of bidirectional 
solitons on water by using integrable Boussinesq surface- variable equation (^). However 
the integrability of the other two systems (0) and (^) or generally the GBQS @ is not 
known. In this paper, we will show that the GBQS is non-Painleve integrable generally 
except for a special equivalent case of (|J). Though the model is non-Painleve integrable, 
it is still possible and useful to give out some interesting solitary wave solutions. In this 
paper, we are concentrate on to study the exact solitary wave solutions of the GBQS 
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under the physically significant boundary conditions 

d n u 



0, 



dx n 

*±oo 



-0, n = 0, 1, 2, ••• . (6) 

x— >±oo 



Because the integrability and the soliton solutions of the model (f|) have been known 
in literature, we always assume ci ^ except for the special cases which will be 
particularly pointed out if it is necessary. 

In the section 2 of this paper, the non-Painleve integrability of the GBQS is stud- 
ied by using the standard Weiss- Tabor-Carnevale (WTC) approach |7|]. A special ex- 
plicit solitary wave solution with a specific velocity selection is given by the truncated 
Painleve expansion. In Sec. 3, after reflecting the problem to find the possible soli- 
tary waves to the possible motions of the Newtonian type quasi-particles moving in 
some proper potential fields, the velocity prohibition and selection phenomena for the 
solitary waves of the GBQS are discussed. For the general solitary wave with allowed 
velocities, an implicit form of the exact solitary wave solutions is given. Some special 
solitary wave solutions are plotted also according to the implicit expression. The last 
section is a short summary and discussions. 

2. Non-Painleve integrability of the GBQS and its explicit exact 

SOLITARY WAVES WITH THE FIRST TYPE OF SPECIAL VELOCITY SELECTION 

2.1 Non-Painleve integrability of the GBQS with c\ ^ 

In the modern soliton theory, the study of the Painleve property @ plays a very 
important role because it can be used not only to isolated out (Painleve) integrable 
models but also to find many other integrable properties such as the Backlund trans- 
formations, Lax pair, Schwarzian form etc@, f|. Furthermore, even if the model is 
non-Painleve integrable, the method can still be used to find some useful things such 
as the special exact explicit solutions |§. Because it has been known that the model 
(fD (i.e., (|5|) with c\ = 0) is integrable and the general real physical model may require 
c\ 7^ 0, say, (||]) and (|3]), we only check the Painleve property of the model with c 7^ 
in this subsection and find a special solitary wave solution from the next subsection. 

The Painleve property of a partial differential equation system is defined as all the 
solutions of the system are free of the essential and branch movable singularities around 
an arbitrary (both characteristic and non-characteristic) singular manifold [|IU|| . 

According to the above definition of the Painleve property, the general solutions of 

the GBQS should have the following expansion around the arbitrary singular manifold 
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(f) = <j){x,y,t) = 



u 



oo 

3=0 3=0 

where {a, (3} should be negative finite integers. Furthermore, if the model possesses the 
Painleve property, there should be at least one primary branch (one possible selection 
of {a, (3, Mo, t>o}) such that enough arbitrary functions (six for the GBQS (0), <fi an d 
five of {uj, Vj}) can be included in the Painleve expansions. 

By using the standard leading order analysis, i.e. substituting u = uo<p a , v = v$$P 
into the GBQS and balancing the leading terms, we can find that there is only one 
possible selection 

a = -2, 13 = -2, u = I2d<j) x <f) u v = 6c 2 (f) 2 x (8) 

for c 2 ^ 0. For c 2 = 0, there is no possible negative integer selection for a and (3. 

Substituting the above expansion(^) with (|8|) to Eq.(^) and vanishing all the coef- 
ficients of the powers (fi k , we can find the following recursion relation for determining 
the expansion functions uj and Vj 

(i-4)(i-6)(l + iK = F y , (9) 
<t> x (j - A)c 2 j{j - 5K- - 120t(j - 4) Cl ^ = F 2j , (10) 

where Fij and F 2 j are all only functions of uq, ui, Uj-i, vo, v±, Vj-i and 
the derivatives of (p. From the recursion relations @ and ([10|), we know that all the 
functions Uj and Vj are fixed by the recursion relations except for those related to the 
resonance points determined by vanishing the coefficient determinant of ([]) and ([TO]): 



(j-4)(j-6)(l + j) 
M-4)c2j(j-5) -120 t (j'-4) Cl 



-120 tCl (i-4) 2 (i + l)(j-6) = O. (11) 

From the expression ([TTD and by checking the resonance conditions at the resonant 
points j = 4, 4 and 6, we know that there are two facts which destroy the Painleve 
integrability of the GBQS with c\ ^ 0. The first one is that there is no primary 
branch at all because the both equations of the QBQS (0) are three order, so the 
primary branch should have six resonant points such that six arbitrary functions can be 
entered into the general Painleve expansion solution. The second one is the resonance 
conditions at j = 4 and j = 6 are not satisfied. The lack of the primary branch of 
the model means some kinds of the logarithmic branch and/or essential singularities 
around the singular manifold <p = may appear for the GBQS (El) with c% ^ 0. 



Actually, substituting u = u (f) a , v = Vq4>P into the GBQS (|5]), we can find that there 
exists a possible non-completely negative selection for the constants a and (3: 

a = -2, /3 = 0, (12) 

for both c 2 ^ and c 2 = cases! The allowance of the selection /3 = means that the 
logarithmic term In <ft should be included in the expansion of the function u to balance 
the leading order terms of the first equation of @. 

According to the above discussions, we can conclude that the GBQS with c\ ^ has 
no Painleve property. In other words, the model is non-Painleve integrable. Though 
the model is non-Painleve integrable, it is still possible to find some useful information 
from the Painleve analysis. In the next subsections we are specially interested to find 
the possible solitary wave solution from the truncated Painleve expansion. 

2.2 Truncated Painleve expansion and explicit exact solitary waves with 
special velocity selections 

According to the discussion of the last subsection, we know that the standard trun- 
cated Painleve expansion has the form 

U = U (f)~ 2 + Uit/T 1 +U 2 , V = V (f)~ 2 + ViCf)' 1 + V 2 , (13) 

where {112, v 2 } is a seed solution of the GBQS (|5|) and usually is taken as constant 
solution to get the single travelling solitary wave solution of the model. 

As usual, substituting the truncated expansion ( |TB| ) with constants u 2 and v 2 into 
the GBQS (|5|) and vanishing the coefficients for different orders of <p, we have 

u = 12ci0 x t , v = 6c 2 (fil, (14) 
12 f 4>x4> tt , 4>t4>xx 



Ul = ~T \JT + + 3< ^ ' ' (15) 



2 A x ^tt tyxh 



^H 2r r-^f- i3 H' (16) 

for the functions «o, vq, u±, v± while the function should be determined from the 
following over-determined system 

u 2 u lx + v lx + uu - ciu lxxt = 0, (17) 

v 2 u lx + v u + u 2 v x + u x - c 2 u lxxx = 0, (18) 

U 0t - Ui<fr t + Ci(2<fi x Ul xt + U lt (f)xx + (Ul<p t )xx - U 0xxt ) 

+U!U lx - U 2 Ui(j) x + U 2 U 0x - V!(f) x + Vox = 0, (19) 
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(u ui) x + 2ci[(p xx u 0t + {4> t u Q ) xx - (j> xx <f> t ui - u lt (fi x ] 

-2u (f) t + [24ci(w t - ui4> t ) x - 2v -u\- 2u 2 u Q )](fi x = 0, (20) 



c 2 [3((f) x ui x ) x + u X (j) 

xxx ^Oxxx ] ~ (U\ + U\V 2 + ViU 2 )(j) x 
(1 + V 2 )u 0x + V 0t + V 0x U 2 + - V X (j)t = 0, (21) 

and 

(m ^i + VqUi) x - 2v (fi t - 2(v 2 u + v u 2 + vtm + u )4> x 

+2[u (p xxx - ^ x {ux^ x ) x + 3(<fr x u ) x ]c 2 = 0. (22) 

To find all the possible exact solutions of the over-determined system is quite difficult. 
However, it is a quite easy work to find the travelling solitary wave from the system 
(|l4])-(p2"|). For the travelling wave solution cfi = (fi(x — ct) = 0(r), we have 

4> x = (fi T , (fi t = ~C(fi T , (23) 

with c being an arbitrary velocity parameter. 

Under the travelling wave solution condition (^), (|l4l)-(|T6|) are simplified to 

M = — 12cic0^, v = Qc 2 (fil, ui = 12cci(fi TT) vi = -Qc 2 (fi TT . (24) 

Substituting (0) into (|T7p and (|T|) yields 

2c 2 cl(fi rTTTT + (2u 2 cci — 2c 2 c\ — c 2 )(fi TTT = 0, (25) 
2ccic 2 (fi TTTTT + [2cci(l + v 2 ) + c 2 (c - u 2 )](fi TTT = 0. (26) 

The equation system ( p5[) and (|26|) is consistent only for 

= c 2 + 2c 2 c 2 (c 2 k 2 -2) c 2 + 2c 2 c 2 (c 2 k 2 - 2) 

V2 4c 2 c? ' U2 2cc 1 c-2 + 2c 1 c 2 (l -ciP) [ ' 

where the constant k is introduced for convenience later. Using the constant relation 
(|27|), the general solution of fl25|) (and fl26|)) reads 

(fi = b + b\T + b 2 r 2 + ai exp(/rr) + a 2 exp(— kr) (28) 

with feo ; &i , b 2 , ao and ai being arbitrary constants. Substituting (|24"D with (|28"D into 
(p~9|) , we find that the constants appeared in (|28| ) have to be fixed by 

b l = b 2 = a x a 2 = 0. (29) 

Furthermore, without loss of generality, the constant a 2 can be taken as zero because k 

is arbitrary (can be taken as both positive and negative). After finishing some simple 
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direct calculation, one can find that (p8|) with ( p9|) solves all other remained equations 

-(H). 

into (|i~3"|), we get a special solitary 



Substituting 

wave 

u 



(|27|), with (0), (g9|) and a 2 
c 2 2 + 2c 2 c 2 (c 2 k 2 -2) 



2ccic-2 + 2cic 2 (l - 
c 2 + 2c 2 c 2 (c 2 A; 2 - 2) 



Cl A; 2 ) 
3 



+ 3cci/c 2 sech 2 

k , 



k 

— [t 
2 K 



TO) 



c 2 k sech 



to 



(30) 
(31) 



4c 2 c 2 2' 

where To = ln(&o) — ln(ai). In general, the boundary conditions expressed in (§) can 
not be satisfied for the solitary wave solution expressed by (p0|) and (|31"D. In order to 
satisfy the boundary conditions (||) for (|^) and (|3lD, we have to fix the constants k 
and c as 



fc 2 



c 2 + 2 



2c 2 ' (2-c 2 ) 
and then the solitary wave solution becomes 

3(c 2 + 2) 



2\ ' °0 



C2 
Cl' 



2^2 



=sech 



2)sech 2 



>/2 Cl (^ + 2) 
4cic 

V / 2c 1 (c 2 + 2) 
4cic 




(32) 

(33) 
(34) 



From the above discussions, the zero boundary travelling solitary wave solution ( p3|) 
(and (|3"4"D) obtained via the truncated Painleve expansion is valid only for the GBQS 
(0) with the condition Cq < 2 (i.e., c 2 < 2ci) and c 2 7^ for c\ 7^ 0. The solitary wave 



solution (H) (and ©) possesses only two special isolated velocities 



Is there 



any other travelling solitary wave with other velocities for the GBQS (|5|)? The answer 
is positive because in the truncated Painleve expansion approach we require all the 
coefficients of the different powers of being zero. This strong condition leads to the 
loss of the generality. In the next section, we discuss this problem generally by casting 
the problem to the possible motions of a Newtonian classical quasi-particle in some 
possible potentials. 



3. Implicit travelling solitary wave solutions and special velocity 

selections and prohibitions 

From the last section, we know that by means of the truncated Painleve expansion 
approach, we can only obtain a special exact solitary wave solution with the boundary 
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condition (|6]). Actually, the travelling solitary waves with different velocities do exist 
for the GBQS (||). Though we can not explicitly write down the exact solitary wave 
solutions with the boundary condition (|6|) for the different velocities, we can find an 
uniform and implicit formula for the solitary wave solutions except for some isolated 
special critical velocities. 

For a travelling wave solution, u = u(x — ct) = u(t), v = v(t), the motion equation 
system (|5|) becomes an ordinary differential equation system 

1 2 

- cu + -u + v = -c x cu TT ^ 

— CV + (1 + v)u = C 2 M rT , 

where the both equations have been integrated once with respect to r and the integra- 
tion constants have been fixed as zero because of the boundary conditions given in (|6|). 
From the first equation of (^), we know that the travelling wave of the field v can be 
simply expressed by 

v = —~{—2cu + u 2 + 2cicu TT ). (36) 

Substituting it into the second equation of (^) we have 

2(1 - c 2 )u - u 3 + 3cu 2 
2ci(cq — c 2 + cu) 

The first integration of ( pTj) reads 



u rr = — TTTT^ 3 , — ■ ( 37 ) 



3 „,2 / „. / r A 



U U I Cn \ U C, 



3ciC ' 2ci V c 2 J C\C \c 2 
-2V{u) (38) 



for c and 



u 4 u 2 



u\ = + — = -2V0(u), t = x, (c 2 ^ 0) (39) 

4c 2 c 2 

for c = 0. To get the relations (p8|) and (|39D , the integration constants have been fixed 
appropriately such that the boundary conditions (|S|) may be satisfied. 

Now, from the expressions (|37|)-(|39D, it is known that to find the possible travelling 
solitary wave solutions of the GBQS (H) with boundary condition (|6|) is equivalent to 
find the possible special motions of a classical quasi-particle moving in the potential 
fields V(u) and/or V0(u) related to the maximum point at u = with the "space" 



variable u and the "time" variable rill 



If the existence problem of the solitary wave is solved, except for some special critical 
cases (see later), the travelling solitary waves with the boundary condition (||]) can be 
uniformly expressed implicitly by 

7b = ±f-riUf, (40) 
' ' -2V(u) 



where V(u) is defined in ( |38"D and the integration constant tq is related to the location 
of the solitary wave. 

To see the possible solitary wave solutions qualitatively, the possible motions of 
the quasi-particle in the potentials V(u) and/or V0(u), it is convenient to study the 
structures of the potentials V{u) and V0(u) at the same time. 

To study the structures of the potentials V(u), we firstly isolate out some special 
and/or critical cases. 

The first special case is same as that in (|32|) . When the velocity c of the solitary 
wave is fixed by (|32"D, the logarithmic term of (^) vanishes and the function u can be 
explicitly integrated out. The result is reasonably same as that obtained in the last 
section via the truncated Painleve analysis. The potential structure at this special case 
is plotted in Fig. la for the special model (0), i. e. c\ — |, c 2 = | and then c = y| 

and c = <\p^- The corresponding soliton solution for the field u is plotted in Fig. lb. 
The second special critical case is that if the velocity is determined by 

c = ±c , (41) 

then the logarithmic term of fl38|) is also vanished. However, in this case there is no 
solitary wave solution with the boundary condition (|6|) except for Co = 1, i.e., C\ = c 2 . 
In other words, the solitary wave solutions with the velocities ±co are prohibited for 
the GBQS (|5|) with c 2 ^ c±, c\ ^ 0. This velocity prohibition property can be observed 
more clearly from ([37]). After substituting (f|l]) into fl3?|) , one can directly see that it 
is impossible to satisfy the boundary condition (|6|) except for cq = 1. While if Cq = 1, 
pip is same as (|3^) and then the related solitary wave solution is also given by (|33|). 
The related potential structure of the potential V(u) with c\ — |, c 2 = | and then 
c = c = </| is plotted in Fig. 2. From Fig. 2, we can see also that when c = c , 
u = is not an extremum of the potential while a solitary wave is corresponding to a 
possible special motion of the quasi-particle related to a maximum of the potential. 
The third special critical case is related to the static solitary wave solution. For the 

static solitary wave solution of the GBQS, two subcases, c 2 > and c 2 < 0, should be 

9 
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Fig. 1b 



1.f> - 




/ 2 - 












/ ^ 




/ 0.5: 





-1 



1tau 



Figure 1. (a). The potential structures for the special velocity selection 



with C\ = \ and Co = 



C2 



(b) The bell shape solitary 



wave solution related to (a). All the figures of this paper have no unit 
because the model is dimensionless. tau = r in all figures also. 
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Figure 2. The potential structures for the critical velocity c = Cq with 
C\ — \ and Co = 



clarified. For the GBQS (||) with c = 0, C2 7^ 0, fl38|) is not valid and has to be changed 

as fl3~9|). The corresponding structure for the potential V0(u) is plotted in Fig. 3a for 

C2 = I > and the corresponding solitary wave solution is plotted in Fig. 3b. 
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Fig. 3a 
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Figure 3. (a) The potential structures for the static velocity c = with 
c 2 = | > 0;. (b) The related static bell shape solitary wave (up solid 
line) and ring shape wave (the lower dashed line). 



From (|39D and (|3~6D , it is straightforward to get that the static solitary wave of the 
GBQS with c 2 > has the form 



u 



±2sech- 



x 



-sech 



x 



(42) 



Fig. 4 shows the structure for the potential V0(u) with c 2 = — | < which includes 



the integrable case (f|). From fl39|) and Fig. 4, it is also easy to see that whence c 2 < 0, 
u = is a minimum of the potential V0(u). So there is no static solitary wave of the 
GBQS with c 2 < and zero boundary conditions because a solitary wave corresponding 
to the motion of the classical quasi-particle in the potential field related to a maximum. 



Furthermore, from (35) with c = and c 2 = 0, we can also know that there is no static 
solitary wave at all. In summary, the static solitary wave with zero boundary condition 
is prohibited for the GBQS with c 2 < 0. Especially, there is no static solitary wave 
with zero boundary conditions for the systems (0) and @ even if @ is integrable. 
The fourth type of velocity selection is related to 

c = ±l, c ^l. (43) 
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Figure 4. The potential structures for the static velocity c = with 



C-2 



< 0. 



Usually, the solitary waves are decay exponentially apart from the solitary wave center. 
However, for the GBQS with C2 7^ ci, if the solitary wave possesses the velocities ±1 
then it decays only algebraically This velocity selection property can be seen from the 
second order differentiation with respect to u of the potential V(u) at the point u — 0. 
The second order differentiation of V(u) for c 2 = 1 reads 



d 2 V(u) 
du 2 



u[Q(c — u — cCq) + 2cu 2 + 3-uc 2 , 
2ci(l -cu-cl) 2 



1). 



(44) 



It is known that the boundary values (u = in our case) of exponentially decayed 
solitary wave solutions are linked with the maxima of the potential However, 
from (^) we know that when the velocity c = ±1, u — is only an inflexion point 
of the potential V(u). It is also known that an inflexion point may be related to a 
rational (or algebraic) solitary wave solution [[TI^]. This type of the algebraic solitary 
wave is similar to some what of the solitary wave at the critical point where the phase 
transition occur for some types of quantum fields and condense matter systems JT^ . 



To see the algebraic decay property of the solitary waves at this critical case, we 



take c = 1, 
becomes 

V{u) 



Cl 



\ and Co = 50. Under this parameter selection, the related potential 



/ u 



15624992502 In I 

V2499 J 3 



u 3 + 1251m 2 - 6252498m = Vl(u) 



139 



-u 



-u = V2(u) 



(45) 

. _ , (46) 
2499 1387778 v ; 

The corresponding structure for the exact potential Vl(u) and the approximate po- 
tential V2(u) are plotted at the same figure, Fig. 5a. From Fig. 5a, it is seen that 

two lines are almost coincide with each other at the plotted region which related to the 
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Figure 5. (a). The potential structures for the potentials Vl(u) and 
V2(u) at the critical velocity c = 1 with c\ = | and cq = 50. (b). The 
corresponding approximately algebraic decayed solitary wave. 



quasi-particle move in. Actually, the approximation (|6|) is quite well up to u ~ ±1000 
while the solitary wave is located only at the region u = ~ 4. So to see the solitary 
wave structure in this special case we can safely use the approximate potential V2(u) 
to calculate the related solitary wave. The result reads 

4998 

W c =l, Cl =0.5, c 2 =100 ~ r 2 _|_ 1251' ' ' 

Fig. 5b. shows the structure of the algebraic solitary wave expressed by (f47|). 

From the analysis of the existence problem of the solitary waves at the critical ve- 
locities, we know that for different velocities the solitary waves may possesses different 
shapes. This type phenomena occurs also in integrable cases. For instance, for the bidi- 
rectional Kaup-Kupershmidt equation, its right moving solitons possess single-humped 
shape while its left moving solitons possess double- humped shape In integrable 
case, the velocity prohibition phenomena are also common. For instance, for the sys- 
tem (^) the static solitary wave with zero boundary condition is prohibited and for the 
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KdV equation ([I]), all the left moving and static soliton solutions with zero boundary 
conditions are prohibited. 

Now the important problem is in addition to the selections and prohibitions at the 
critical points whether there exist further prohibitions for some other velocity regions. 
To solve this problem, we restrict ourselves for c\ > (because most of the real physical 
systems, say, @ and (|3|), possess this property) and c > (because the symmetry 
property {c, u] <-> {— c, — u] of ([57])). 

To find the existence conditions of the solitary waves of the GBQS (|5]) with zero 
boundary conditions (f)|) for noncritical cases, we can use the maximum condition, 
< and the real condition of the potential V(u). 



d 2 V(u) 
du 2 



u=0 

The maximum condition of the potential V(u) at u = reads 



du 2 

and the real condition of V(u), reads 



,, - T "" -4^T<0, (48) 
u=0 ci(c 2 -c^) 



cu , , 

1 > 0. (49) 



c 2 



After finishing the detail analysis with help of the conditions (fig) and (49), we find 
the following six different structures for the potential V(u) for c\ > 0: 
Case 1. 

< c 2 < c 2 < 1. (50) 

When c 2 , the square of the velocity parameter, is located in the range (0, cq) with 
c < 1, the related potential V{u) possesses the structure as shown in Fig. 6a for 
c > 0. From Fig. 6a we know that there exist a special motion for a classical quasi- 
particle moving in this potential related to the maximum at u = 0: At the beginning 
(r = — oo), the quasi-particle is located at the peak center u = 0, as time "r" increases, 
the quasi-particle "roll" down the hill up to the um = u am point and then return back 
to the original point u — at "time" r = +oo. This type of special motion of the 
quasi-particle is just related to the solitary wave solution of the GBQS @ with ( |50|) . 
From Fig. 6a we know also that the quasi-particle can roll only to right (the positive u 
direction) that means the right moving solitary wave (c > 0) possesses bell shape for the 
field u. By using the invariant transformation of (|37|), {u, c} — > {—u, — c}, we know 
that the left moving solitary wave (c < 0) possesses ring shape. The corresponding 
exact solitary wave with the same parameters as shown in Fig. 6a is plotted in Fig. 
6b. 

14 



Fig. 6a 
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Figure 6. (a) The typical structure for the potential of the quasi- 
particle located in with c\ > and c>0,0<c 2 <Cq<1 and 
{c = c = 0.4, ci = |}. (b) The corresponding solitary wave 
solution related to (a). 

Case 2. 

cl<c 2 < 1. (51) 

In this case, the related potential V(u) possesses the structure as shown in Fig. 7 for 
c > 0. From Fig. 7 we know that there is no solitary wave solution with zero boundary 
condition in this case because u = is only a minimum of the potential. If the quasi- 
particle stay at u — at the beginning then it can only be stayed there forever. In 
other words, the solitary waves of the GBQS ([D with the velocities Cq < c 2 < 1 and 
zero boundary condition @ are totally prohibited. 
Case 3. 

eg < 1, c 2 > 1. (52) 
The structure of the potential V(u) for c 2 > 1 with Cq < 1 possesses the form as shown 
in Fig. 8a. Similar to the case 1, there exist the right moving bell shape solitary waves 
and the left moving ring shape solitary waves for the field u. The right moving bell 
shape solitary wave solution related to Fig. 8a is plotted in Fig. 8b. 
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Fig. 8a 




Fig. 8b 
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Figure 8. (a). The typical structures for potential V(u) with c\ > 0, 
c > 0, c\ < 1, c 2 > 1 and {c = 4=, c = 2, c x = \}. (b). The 
corresponding solitary wave related to (a). 



In the first three cases, cq is less than 1 and the special models (0) and @ satisfy this 
condition. When this condition is not satisfied, we have three other different potential 
structures. 
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Fig. 9b 
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Figure 9. (a). The typical structure for potential of the quasi-particle 
located in with c\ > 0, c > 0, Cq > 1, c 2 < 1 and {c = 2, c = 0.6, C\ = 
|}. (b). The corresponding solitary wave related to (a). 



Case 4- 

< c 2 < 1 < eg. (53) 

The corresponding potential structure of the fourth case (|53|) is plotted in Fig. 9a. 
Different from the first and third cases, from Fig. 9a we know the quasi-particle located 
the u = hill may roll down the hill to both the left and right sides and finally return 
back to the hill. So in this case, the right moving solitary waves may have both bell and 
ring shapes and the left moving solitary waves possess the same property. The right 
moving bell shape solitary wave solution related to this case with the same parameters 
as Fig. 9a is plotted in Fig. 9b. From Fig. 9a, we can know also that for the right 
moving solitary waves, the bell shape solitary waves possess larger amplitudes and for 
the left moving solitary waves, the ring shape solitary waves possess larger amplitudes. 
Case 5. 



1 < C 2 < Cq. 
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(54) 



Figure 10. The typical structures for potential V(u) with c\ > 0, c > 0, 
1 < c 2 < Cg and {c = 2, c = 1.6, ci = -}. 
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Figure 11. The typical structure for potential V(u) with c\ > 0, c > 0, 
1 < Cq < c 2 and {c = 2, c = 2.1, c x = |}. 



When c 2 is located in the range (1, cq), the corresponding potential structure possesses 
the form as shown in Fig. 10. From Fig. 10 we know that u = is only a minimum of 
the potential. So similar to the discussion of the case 2, we can conclude that there is 
no zero boundary solitary wave for the GBQS with the condition (54). 
Case 6. 

1 < c 2 < c 2 . (55) 

For the last case (|55|), the related potential structure is shown in Fig. 11. From Fig. 11 
we know that u = is really a maximum of the potential, the quasi-particle may roll 
down the u = hill in both sides, however, it can not be back! That means the fast 
moving (c 2 > Cq) solitary waves for the GBQS with the condition Cq > 1 are completely 
prohibited also. 
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Combining the critical cases, we summarize the conclusions for the possible solitary 
waves of the GBQS (f|) with zero boundary conditions (f|) and c\ > in the following 
table. 
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Table 1. Existent and prohibited regions of the solitary waves. 



model 


velocity 


solitary wave 


shape 


C 2 < Ci 


c< -1 


yes 


ring shape {u < 0) 


C 2 < Ci 


— 1 < c < — Co 


no 




< c 2 < ci 


- c < c < 


yes 


ring shape 


c 2 < 


- c < c < 


no 




< c 2 < ci 


c = 


yes 


both ring shape and bell shape 


c 2 < 


c = 


no 




< c 2 < ci 


< c < c 


yes 


bell shape (u > 0) 


c 2 < 


< c < c 


no 




C 2 < Ci 


Co < c < 1 


no 




C 2 < Ci 


c > 1 


yes 


bell shape 


Ci < c 2 


c< -1 


no 




Ci < c 2 


c= -1 


yes 


algebraic, ring shape 


Cl < c 2 


-1 < c< 1 


yes 


both ring shape and bell shape 


Ci < c 2 


c = 1 


yes 


algebraic, bell shape 


Cl < c 2 


c> -1 


no 





4. Summary and discussions 

In summary, using the standard WTC's Painleve PDE test, we prove that the GBQS 
is non-Painleve integrable except for the special case for ci = 0, i.e., (^). For the non- 
Painleve integrable GBQS with zero boundary conditions @, the truncated Painleve 
expansion approach leads to only a special sech 2 shape solitary wave solution with 
a special velocity selection. To find all the possible travelling solitary wave solutions 
of the GBQS with ci > for all the possible model parameter regions, we map the 
problem to find the possible motions of a Newtonian classical quasi-particle moving 
in some possible potentials. After considering all the possible motions of the classical 
quasi-particle in the possible potentials related to the maxima of the potential at u — 0, 
all the possible travelling solitary wave solutions related to zero boundary conditions 
are found for all the possible velocities and the model parameter ranges with ci > 0. 

Similar to the integrable cases such as the KdV equation, the solitary waves at 
some special ranges, the travelling solitary waves with zero boundary conditions are 
completely prohibited. For the KdV equation, all the left moving and static solitons 
with zero boundary conditions are prohibited. For the GBQS with ci > 0, there are 
three different cases. For the first case, < c 2 < c\, both the faster and slower moving 
solitary waves are allowed while the solitary waves moving in the "middle" velocities 
^ = Cq < c 2 < 1 are completely prohibited. For the second case, < — c 2 < c±, only 
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the faster moving solitary waves are allowed while all the slower moving solitary waves 
c 2 < 1 are prohibited. For the third case < c\ < C2, contrary to the second case, 
all the slower moving solitary waves are allowed while all the faster moving c 2 > 1 
solitary waves are completely prohibited. For the first two cases the solitary waves at 
the critical velocities c = ±1 are prohibited while for the third case, there are algebraic 
solitary waves at the critical velocities c = ±1. For the first two types of GBQS, 
the right moving solitary waves possess bell shape and the left moving solitary waves 
possess ring shape. For the third type of GBQS, the right moving bell shape solitary 
wave possess larger amplitudes than those of the right moving ring shape solitary waves 
(with the same value of the velocity) while for the left moving solitary waves, the ring 
shape solitary waves possess larger amplitude. 
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